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£Nj . In this series of three papers, we generalize the derivation of dual photons and monopoles by 

Polyakov, and Banks, Myerson and Kogut, to obtain approximative models of SU(2) lattice gauge 
theory. Our approach is based on stationary phase approximations. 

In this second article, we start from the spin foam representation of 3-dimensional SU(2) lattice 
gauge theory. By extending an earlier work of Diakonov and Petrov, we approximate the expectation 
value of a Wilson loop by a path integral over a dual gluon field and monopole-like degrees of 
freedom. The action contains the tree-level Coulomb interaction and a nonlinear coupling between 
dual gluons, monopoles and current. 
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I. INTRODUCTION 
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It is well-known that in the strong-couplingregime of pure SU(N) gauge theory, confinement can be derived from an 
expansion in strong-coupling graphs [3, ^ arises from three steps: an expansion in characters, an integration 
over the group variables, and an expansion in the inverse coupling. It is less known that the first two steps can be also 
»£m done explicitly without any strong-coupling expansion: then, one obtains a non-perturbative sum over graphs that is 
I 1 exactly equivalent to the original lattice gauge theory 0, In this context, it would be misleading to speak of 
Oh, "strong-coupling" graphs, so we follow the paper Q and refer to them as spin foams 1 . 

Physically, the spin foams may be viewed as worldsheets of flux lines. Thus, the spin foam representation is an 
approach to gauge theory where everything is formulated in terms of flux lines instead of fields. 

Up until recently, the strong-coupling expansion was the only example of an analytic computation with spin foams 
of gauge theory. Due to the strong coupling, fluctuations of the spin foam surfaces are strongly suppressed. The 
dominant contribution to the Wilson loop comes from the minimal surface that is spanned by the loop, and this gives 
the area law. The problem about this argument is that it requires a size of the lattice spacing that is comparable 
to the length scale of confinement. Therefore, it cannot serve as a completely satisfactory explanation. It would be 
preferable to have a derivation in the continuum limit that can produce both the short- and large-distance potential in 
a coherent fashion: i.e. with the Coulomb potential dominating at short distances, and the confining potential taking 
over at a certain length scale above the lattice cutoff. 

This naturally suggests the following question: how can we perform an analytic computation with spin foams when 
the lattice spacing is small? In the continuum limit, the spin foam surfaces undergo complicated fluctuations and it 
is not clear how one could sum over them. 

A main motivation for the string approach to gauge theory comes from a very similar question (for a review, see 
e.g. 0, [l(| EH now can one nn d a continuum representation of gauge theory that implements the surface- like 

(or stringy) aspects of the strong-coupling expansion? 

The answer to our questions is known for the analogous problem of U(l) lattice gauge theory: there the counterpart 
of the spin foam formulation could be called a "charge foam" representation. It is given by a sum over branched 
surfaces that are labelled by U(l) representations (i.e. charges) and constrained by the Gauss law. In this context, 
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1 This term was originally coined in the quantum gravity literature Q. 
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the problem of confinement is solved by going to yet another representation: as was shown by Banks, Myerson and 
Kogut, the charge foam formulation can be transformed exactly to another representation that has dual photons and 
monopoles as its degrees of freedom. It was derived earlier by a different method by Polyakov [l3|, [3] • The photon- 
monopole representation allows for an analytic computation of the static potential between charges, and provides a 
direct explanation for confined and deconfined phases of the theory. For suitable values of the dimension and coupling, 
the monopoles condense along a string between the charges and thereby create a confining potential between them 

[flElIIEOa. 

The example of U(l) fosters the hope that one could generalize this scheme to non-abelian gauge theories: is there 
an analytic way to deduce analogues of the photon-monopole representation for SU(N) lattice gauge theory? In the 
three papers of this series, we address this question and derive approximative gluon-monopole models from three 
representations of SU(2) lattice Yang-Mills theory: from the BF Yang-Mills representation in dimension 3 and 4 [Isj |. 
from the spin foam representation in d = 3, and from the plaquette representation in d = 3 p^ |. 

In this paper, we will start from the spin foam representation in 3 dimensions and attempt to generalize the 
aforementioned transition from charge foams to dual photons + monopoles. We thereby extend an earlier work by 
Diakonov and Petrov [20l |. where it was already suggested how dual gluons arise from spin foams. 

The logic of the derivation is analogous to that of paper I: every step that was made within the BF Yang-Mills 
representation has a counterpart in the spin foam picture. It is non-trivial, however, to work out the details of this 
correspondence. The resulting model is similar to what we obtain in paper I and III: it contains a tree-level Coulomb 
interaction, and a nonlinear coupling between monopolc-like excitations, dual gluons and current. 

The paper is organized as follows: after setting our conventions for SU(2) lattice gauge theory (sec. |TT|, we will 
review its spin foam representation in section IIIII In section IIVI we explain how notions of BF Yang-Mills theory 
translate into the spin foam representation. The main part is section |Vj where we derive the representation in terms 
of dual gluons and monopole-like excitations. The final section contains a summary and discussion of the results. 

For conventions on lattice quantities, see the introduction of paper I [181 ]. 



II. SU(2) LATTICE YANG-MILLS THEORY IN 3 DIMENSIONS 

The partition of function 3-dimensional SU(2) lattice Yang-Mills theory is defined by a path integral over SU(2)- 
valued link (or edge) variables U e on the lattice k: 



Z = J Ql^e) exp ( - S f( W f) 



(1) 



The face (or plaquette) action Sf depends on the holonomy Wj_ around the face. As in paper I, we choose Sf to be 
the heat kernel action (for more details on the definition, see 2l|). The heat kernel action has a particularly simple 
expansion in terms of characters, namely, 

exp(-5 / (W / )) =^(2i + l)e-f^' +1 ) Xi (W / ). (2) 
j 

The coupling factor j3 is related to the gauge coupling g via 

4 1 

P = — + o- (3) 



III. SPIN FOAM REPRESENTATION 



Partition function 



The SU(2) lattice gauge theory can be rewritten in different representations. One of them is the representation in 
terms of spin foam sums (d, 0, Q; for a review, see e.g. [22j]). It results from the original path integral in two steps: 
first we expand the plaquette actions in gauge- invariant loops (a.k.a. characters), so that 



Z= I ( n dC/e) £ II (2i/ + l)Xi / (W» e-?^+i) ] (4) 
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Figure 1: Even cubes of the lattice re. 

We may view this as a first-order formulation, since it has two sets of variables: the connection variables U e and the 
spin assignments jf . 

In the second step we integrate out the connection variables U e . This can be done exactly and leads to a sum over 
configurations that we call spin foams. Each spin foam consists of an assignment of spins jf and intertwiners I e to 
faces and edges of re respectively. One may think of it as a branched surface that carries spin and intertwiner labels 2 . 

In general, there are several ways of organizing the integration steps, and each of them results in a different 
form of the spin foam sum. Here, we will use an integration scheme that was proposed by Anishetty, Cheluvaraja, 
Sharatchandra and Mathur H: it is specific to 3 dimensions, and leads to 6j-symbols in the amplitude. In this 
case, the intertwiners are encoded by spin assignments j e to edges of K, and the spin foams reduce to assignments of 
spins jf and j e to faces and edges respectively. To describe the amplitudes, it is convenient to use the following dual 
construction: imagine that we divide the cubes of re into sets of "even" and "odd" cubes, forming a "checkerboard" 
as indicated in fig. Fig.[T] Dually, we have a set of even and odd vertices in the dual lattice re*. Take the odd vertices, 
and connect each pair of odd vertices within the same face by a diagonal. In this way, the dual lattice re* turns into 
a triangulation T (see Fig. Faces of n are dual to edges of k*, and the edges of re are dual to the edges we added 
to re* to get T. Thus, a spin foam can be described by assigning a spin j e to each edge e of the triangulation T. The 
spin foam sum takes the following form: 

Z= E (UW' + v) (liA (il (-l) 2 -e-|>^)) (5) 

{je>T VeCT / VtCT / VeCK* / 

Each configuration {jej-T is an assignment of spins j e to edges e of T such that for each triangle of T the spins satisfy 
the triangle inequality. The edges of T belong to two groups: edges that are identified with edges in re*, and diagonal 
edges that were added to re* in order to form the triangulation T. In the amplitude, every edge contributes with 
the dimension 2j e + 1 of the representaiton j e . In addition, edges of re* give a sign factor and an exponential of the 
Casimir. For each tetrahedron t, we get a 6j-symbol 

fr- — je 

A t = { h - /:; 1 

1 J4 % J6 J 




(6) 



where the spins ji, ji and j'3 are read off from any triangle in the tetrahedron. The spins j'4, j'5, and je are the 
spins on the edges opposing those of ji, j'2 and j'3. Consider, for example, a tetrahedron that lies in the center of a 



2 Spin foams are essentially the same as the strong-coupling graphs of the strong-coupling expansion [J. Their appearance, however, is 
not tied to any expansion in the coupling, so we prefer to use the term "spin foam" , which was coined in the quantum gravity literature 
Q. For a review of gravity spin foams, see e.g. |2.'I . 



Figure 2: Tetrahedra of the triangulation T that form a cube of the dual lattice k* 



dual cube: if we label its edges by spins as in fig. Fig. [2j the associated amplitude is precisely expression ([6]). The 
6j-symbol is equal to the value of a tetrahedral spin network (see the right-hand side of eq. ©): the labelling of this 
spin network is dual to that of the tetrahedron in the sense that vertices of t correspond to cycles in the spin network, 
and triangles of t correspond to vertices 3 . 

Wilson loop 

How do we compute expectation values in the spin foam representation? Consider a Wilson loop C in the rep- 
resentation j, where C has no self-intersections. In the original formulation, the expectation value of C is given 
by 

(tvjWc) = J (jj.^e) cxp ( - £ S f (W f )) Xj (W c ) ■ (7) 

Wc denotes the holonomy around the Wilson loop. To describe the associated spin foam sum, we have to modify the 
triangulation T: consider all faces of the dual lattice k* which are dual to edges of the Wilson loop. In each of these 
faces, we add a second diagonal edge. Together with the original diagonals, these edges define a "hexahedron" h in 
the center of the cube. We call the resulting new complex T" (see Fig. [3]). 
Then, the spin foam sum can be specified as follows: 

i^wc) =\ e (n + !)) ( n A ( n a9 A ( n (-^ (8) 

{j e } T , \eCT> J \tCT' ) \hCT' ) \eCK* / 

Let us explain the differences to formula ([5]): each configuration {j e }T' is an assignment of spins e to edges of T' 
such that 1. for each triangle of T' the spins satisfy the triangle inequality, and 2. for every double diagonal the spins 
satisfy the inequality 

\ji-j\<fi<ji+j- (9) 
In addition to 6j-symbols, we now have 9j-symbols — one for each hexahedron h in T 1 where the Wilson loop passes 



So that 6j-symbol and spin network are equal, one has to specify suitable phase conventions for the intertwiners of the spin network. 



Figure 3: Modified triangulation T' with double edges where the Wilson line passes through. 



through 4 : 



4 j 



± 



3l 32 J3 

3b 3 35 
36 32 34, 




(10) 



The 9j-symbol is equal to a spin network with an additional edge of spin j. 

IV. RELATION BETWEEN BF YANG-MILLS THEORY, SPIN FOAM SUM AND GRAVITY 

To derive an approximate model of the full lattice gauge theory, we use a similar argument as in paper I. This 
time, however, we will apply it at the level of the spin foam representation. In paper I, we considered the variation 
of the BF Yang-Mills action on the lattice, and argued that for high enough momenta k, stationary phase points are 
given by configurations where the connection is flat and the B-field satisfies a Gauss constraint. We proceeded as if 
this was valid for all energy scales, and restricted the entire path integral to these stationary phase points: i.e. to flat 
connections and solutions of the Gauss constraint. The latter are expressed by (d — 3)-forms on the dual lattice K* , 
so in d — 3 we arrive at a path integral over a scalar field. 

To repeat this argument, we need to know how the variation of the JSF-term and the £? 2 -term translates into the 
framework of the spin foam representation. In paper I, we rewrote the SU(2) lattice Yang-Mills as a BF Yang-Mills 
theory on the lattice: 




(11) 



4 For details about sign factors, we refer the reader to Diakonov & Petrov's paper |2Cj . 
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It can be regarded as a deformation of BF theory on the lattice: 



Z BF — 



n 



( n dO / f n d 3 6^ e -i) I 



1^x^1/2 



sin(|w a;AIly |/2) / 

Both of these path integrals have an equivalent representation as a spin foam sum, namely 



{je} VeCT / \ e £ K » / VtCT 

{je} VeCT / \ e ^ K » / VtCT 



n e -» 

\eCK* 



Je(je + 1) 



We see that the difference lies in the deformation factor 



exp ( --gletie + K 



(12) 

(13) 
(14) 

(15) 



which is the spin foam analogue of the deformation term B 2 in (jlljl . The BF-term, on the other hand, should be 
encoded in the tetrahedral amplitudes appearing in (|13p and (| 14|) . In order to make this correspondence concrete, we 
need to observe two things: 1. 3d SU(2) BF theory is equivalent to 3d first-order gravity. 2. In the large spin limit, 
the tetrahedral amplitudes are related to the action of 3d gravity. 
In the continuum, the BF action is given by 



S- 2 |d 3 x e ptlv B xp -F x , v . 



(16) 



The Lie algebras of SU(2) and SO (3) are isomorphic under the map (of basis elements) 

su(2)3icr a /2 i ► J a eso(3), a = 1,2, 3, (17) 

where 

y a) b c ._ e ab c (18) 

Thus, we can interpret the forms A a and F a also as components of so(3)-valued forms A and F. Then, we get 

F b c = F a (J a ) b c = F a e a b c , (19) 

or 



F a = - e a bc F b 
2 



(20) 



By inserting this in P^|) . and identifying B a with the co-triad e a , we see that the BF term is the same as the 
Hilbert-Palatini action of 3d gravity: 



S = ~ I d 3 x ie^ v e nhr _e a n F hc 



^abc & p,v 

Let us now consider the tetrahedral amplitudes in (fT3|) : they are given by 6j-symbols 

_ ,J Jei Je 2 Je 3 
Je 4 Je 5 Je 6 



(21) 



(22) 



Ponzano and Regge [2J| interpreted the 6 spin assignments j e of a tetrahedron t as lengths j e + 1/2 of its edges e, 
and noticed that for non-degenerate tetrahedra and large spins the amplitudes 



A t = 



Jei Je2 Je3 
Je 4 ]es 3e& 



(23) 
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have the asymptotics 



+ e 



-iRt 



)• 



(24) 



Here, V is the volume of the tetrahedron, and Rt is the contribution of the tetrahedron to the Regge action — a discrete 
version of / d 3 Xy/g R for piecewise fiat geometries [25j : 



tCT 



tCT 



£C? o + l/2)0 te + 7r/4 



,eCt 



(25) 



#t e stands for the dihedral angle at the edge e C t in the tetrahedron t. This motivated Ponzano and Regge to propose 
the sum (1141) as a path integral quantization of 3-dimensional Euclidean gravity: the spin configurations describe 3d 
geometries and the 6j-symbols provide the weighting with the gravity action. 

In this way, we see how the f?F-term in (fTTj) and the tetrahedral amplitudes in (fT"3|) are related: both the BF-term 
and the tetrahedral amplitudes implement a weighting with a gravitational action. In the BF Yang-Mills path integral, 
we have the action of first-order Hilbert-Palatini gravity where the connection is an independent variable, while in 
the spin foam sum the connection is integrated out and the gravity term is that of second-order Einstein gravity. 



V. DUAL GLUONS AND MONOPOLE-LIKE EXCITATIONS FROM SPIN FOAMS 



In this section, we start from the spin foam representation of the Wilson loop average, and derive an approximate 
representation in terms of a dual gluon field and monopole-likc variables. We thereby extend an earlier argument by 
Diakonov & Petrov [2fJ, where it was already suggested how dual gluons arise from spin foams. The line of reasoning 
stands in analogy to what we did in paper I within the BF Yang-Mills representation. Technically, however, the 
present analysis is more complicated and we are forced to make various heuristic assumptions. The constraints, for 
example, give more solutions than in paper I — a feature that could be important for color screening. 

The derivation proceeds in four steps: 

1. First we will apply the Poisson summation formula. We thereby trade the discreteness of spins for additional 
degrees of freedom which are similar to monopoles of U(l). 

2. We use the asymptotic formula for 6j symbols. 

3. By applying a stationary phase approximation we obtain a constraint on spin foams: it requires flatness outside 
and torsion along the Wilson loop. 

4. Solving the constraint yields the dual gluon degrees of freedom. 



Poisson summation formula 



Let us start from the spin foam sums of formula (JSJ) and |(5J). We want to consider the continuum limit, so the 
lattice constraint a is small and (3 is much larger than 1. This means that the dampening due to the Casimir factor 
(|15p is weak, and that typical spins are much larger than 1/2. We therefore approximate 

2j e + l«2j e , j e (j e + 1) » jl , ^1. (26) 

A short calculation shows that we can apply the Poisson summation formula similarly as for U(l), giving us 5 




5 We omit constants that will drop out in expectation values. 
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and 



(trjWcH- f°° (U Jedj e ) E 

x f II A ( II 4 J ) ex P f-i 2 2*i. -E^+E ^ 

VtCT' / \hCT' / \ ecu* eC«* ^ eCT' 



' ' ' ] . (28) 



It is assumed that we replace the discrete functions A t and by suitable continuous functions of spin, and that 
these functions have only support on spins that satisfy the triangle inequalities. 



Stationary phase approximation 

Based on our previous considerations, we identify the 6j and 9j symbols in (f2"5)l as the parts of the amplitude that 
correspond to the BF term and the source current of the Wilson loop. In order to determine the stationary phase 
points of the sum, we proceed in two steps: we will first consider the case without Wilson loop (i.e. the partition 
function |27p) and then discuss how the stationary points are modified, when the Wilson loop is included. Without 
Wilson loop the argument is the same as in Diakonov & Petrov's paper. We repeat it here for completeness. 

If we assume that large spins dominate, the amplitudes of non-degenerate tetrahedra are described by the Ponzano- 
Regge formula By plugging this into (|2"T|) , we obtain 



Z=r( II j^e) EE 
J ° WT J {m e} { s t} 

( n ww exp (' is t R t) ) ex p ( _i E 27r ^ e ~ E ~k je + E 4?ri > me I 



{m c } {s t } 

" " ' " ' (29) 



We sum over signs s t — ±1 for each tetrahedron, due to the two terms of opposite phase in formula 

Let us try to determine the stationary phase points of this expression. If we ignore the polynomial volume factors, 
the relevant exponent is 



*^E 



St 



tcT 



Eie te + 7r/4: 



eCt 



i E 27Tj e - E |-?'e + E 47ri ^ 



= i e je e e + e ^4 - E I ^ 2 + E 4?ri > ( 3 °) 

eCT eCT eCK* " eCT 



where we set j e + 1/2 w j e and defined the angle 



The sum over alternating signs St leads to rapid phase oscillations, so we expect that a stationary phase point can 
be only reached if the signs St are either all +1 or all — 1. In that case, <d e is the deficit angle at the edge e: an edge 
in k* is contained in 4 tetrahedra, and a diagonal edge lies in 6 tetrahedra; hence, for a flat geometry, the sum of 
dihedral angles at e is Att — 2tt = 2tt and 6n — 2tt = At: respectively. 



The first term in (|30[) corresponds to the continuum expression ^/detg R. In the spin foam geometry, lengths are 
given by spins, so the metric is of the order g ~ j 2 . Since the Ricci scalar R goes roughly like g~ 1 g~ 1 ddg, this gives 
us 

jO~ Vdet^i?- yfdrtg~ g^g^ddg ~ (f)^j- 2 j- 2 ddj 2 ~ ddj . (32) 

For Fourier modes j(k) of the spin field, where k is sufficiently high, the phase oscillations in should be dominated 
by the curvature term. Therefore, a stationary phase point is reached when e is zero for all edges, i.e. when the 
spin foam geometry is flat. This is analogous to the argument in paper I, where we concluded that for high momenta 
stationarity is attained by flat connections. Like there, we will proceed as if this was valid for the whole momentum 
range, i.e. we will restrict the entire sum (|27p to flat spin foams. 
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How does the situation change, when we include the Wilson loop? Then, the analysis is less clear, since we have no 
asymptotic formula for the 9j-symbols associated to "hexahedra". We know, however, that the spins ji are typically 
larger than the fixed j of the source, and that the 9j-symbol is zero unless ji and j satisfy the triangle inequality 
(see eq. (TTU10 . The simplest scenario would be that the 9j-symbol has just the effect of ensuring the triangle inequality, 
and can be otherwise approximated by 6j-symbols: that is, if 



the modulus 



< fi < ji+j ■ 



A* 



Jl 


32 


.73 


& 


j 


35 


36 




Ji 



is approximated by the asymptotics of 



and otherwise 





32 


33 1 




35 


k J 






0. 



(33) 



(34) 



(35) 



(36) 



If this is true, what will be the stationary points? Far from the Wilson loop the stationary spin foams should be flat 
according to our previous argument. Directly at the Wilson loop this argument does not apply, since we do not know 
the influence of the phase of AJ . We only know that the additional condition (J33J) has to be satisfied. 

Since we lack a more precise analysis, we will assume that the stationary points are given by all spin foams that 
are flat outside of C and meet condition (|33|) along the Wilson loop. By "flat outside of C" we mean the following 
property: one can 

• remove the union H of all "hexahedra" h along C from T", and 

• map the remaining spin foam into Euclidean R 3 such that 

• the induced metric on T'\H is flat, and 

• spin assignments coincide with edge lengths determined by this metric. 

Such a map sends every edge e of T" into a difference vector b e in R 3 and the length of b e agrees with the spin: 



\be\ = V^ = j e + l/2 



(37) 



Note that for any closed curve around a triangle, the corresponding image vectors b e close as well. If we view the 6 e 's 
as a 2-chain b on T", it satisfies db = for all faces outside the dual *C. Around faces / C T' that are dual to C the 
vectors b e do not close, in general, since b e and b e i may have different lengths (see Fig. They are only restricted 
by the triangle inequality ([33]) : 



3 



< 



\K\ 



< 



\b a \+3 



(38) 



We can regard this failure of vectors to close as a defect in our otherwise flat spin foam geometry: the technical term 
for this deviation is torsion, and the 2-form db measures the amount of torsion in the spin foam. The Wilson loop C 
represents the worldlines of quarks, and that we get nonzero torsion along them fits well with the fact that classically 
a fermion current is a source of torsion (26J. 

The torsion along the loop is constrained by the triangle inequality. We can satisfy it by interpreting the spin j in 
(|38|) as the length of a third vector, and make all three vectors sum up to zero. More precisely, condition (f38|) is met 
if there is a vector 



J = jn, 



n e 



1 



(39) 



such that either 



(I) b e -b e , = J, or (II) b e + b e , = J . 



(40) 
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Figure 4: Torsion vector Jf on faces / dual to C. 



Suppose that the edges e and e' are oriented such that df = e — e'. Then, possibility (I) means that 



db 



f 



J- 



(41) 



and the magnitude of torsion is determined by the spin j of the Wilson loop. The equation looks much like the abelian 
Gauss constraint of U(l) lattice gauge theory. Solution (II), on the other hand, leads to a torsion 



db, 



2b P 



J« 2b e 



(42) 



and the abelian Gauss law is violated. In that case, the torsion would be of the order of the typical spins in the spin 
foams, which is large. In the following, we will drop possibility (II) and only work with the solution (I). There are 
two motivations for this: the first one is simply that we do not know how to analyze condition (II). Secondly, we are 
not sure that the solution (II) would appear in an exact analysis of 9j-symbols. We would expect that the magnitude 
of torsion is determined by the representation of the fermions, and not by the representation of typical fluctuations. 

Even with this assumption the Gauss law (|4ip is not the only way of solving the triangle inequality. The triangle 
inequality states a condition on the length of the vectors b e and b e i, but does not necessarily imply that there is a 
third vector of length j with which they sum up to zero. We can, for example, take b e and b e / to be the same and 
meet the inequality. If we combine this with the condition *d& = = d*b outside the loop, we can construct a solution 
by using a tube-like surface S that ends on the Wilson loop C, and carries some constant vector v € M 3 as "color 
charge" : 



* b = v S 



(43) 



It corresponds to a spin foam with surface S and uniform spin Recall now that the sum over monopole-like 
variables has the effect of restraining spins to be integers or half-integers. As a result, the spin foam surface (|43|) is 
suppressed unless \v\ e Z+/2. Then, however, the triangle inequality 



-J\ 



< 



< 



\v\+j 



(44) 



is only met if the spin j is integerl In other words, the tube-like configuration (|43|) is only allowed when the Wilson 
loop has integer spin. 

Such type of spin foams are precisely the diagrams that are responsible for color screening in the strong-coupling 
regime. We suspect them to be also the cause for color screening in the continuum limit. Our argument shows that 
these spin foams appear as solutions to our flatness and defect condition. We also see that the Gauss law (|4Tj) has a 
smaller set of solutions and cannot capture this non-abelian effect. 

For simplicity, we will not deal with these additional solutions and instead just work with the naive condition (141 1) . 
In that case, the deviation from db = can be described by an R 3 -valued 2-chain J that has length j along *C, and 
is otherwise zero. Thus, we can express the flatness and torsion condition by a single condition on the 2-chain b: 



db = J . 

Since C has no self-intersections and dJ = d 2 b = 0, the vector Jf is always the same along C, i.e. 



(45) 



Jf 



j n , / dual to C . 
, otherwise , 



(46) 
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for some unit norm vector n. This provides the desired characterization of stationary points: a spin foam is stationary 
if there is an R 3 -valued 1-chain b on T" and a current J of the form such that db — J and j e — \b e \ for all edges 
e C T. 

We now constrain the path integral ([28")) to these stationary phase points 6 . The constraint has the general solution 

b = dip + b, (47) 

where ip is an K 3 -valued scalar field on the dual lattice k* and b xll is a particular solution. The latter can be chosen 
as 

b = jn*S (48) 
for any surface S bounded by C. Thus, we can make a change of variables 



(Ujedje) ... - / (ndVWl 

JO V e / JM 3 \v J 



(49) 



and express the constrained path integral by an unconstrained path integral over the scalar field. The Jacobian is 
given by 

!1 , e incoming at v , 
-1, e outgoing at v, (50) 
, otherwise . 

When making this change of variables, we have to observe that the map from ip to j has global degeneracies: a global 
translation and rotation of the values of <p in R 3 does not affect the spins. We fix the translational degeneracy by 
removing the Fourier mode of momentum zero from the (^-integration. The rotational invariance produces a finite 
and field-independent volume factor (the volume of the rotation group) and drops out in expectation values. With 
this, the path integral over <p takes the form 

X ( II j ( II A h ) E 271 ^ - E \ ( d ^+b e ) 2 + 47Ti \dp e + h e I TTlA ■ 

\tcT< ) \hcT' ) V ecu* eCs* 1 eCT' I 

(51) 

We have included an integration over the vector n in b, since its choice is arbitrary. 

The restriction to flat spin foams simplifies the contribution from 6j- and 9j-symbols. Recall that the asymptotic 
expression for 6j-symbols contains volume factors Vt and phase factors exp(iR t ). The latter combined with the sign 
factors (— l) 2jc to give deficit angles 0t in the exponent. In our integral over flat spin foams, these deficit angles 
are all zero. For simplicity, we will also drop all factors with a polynomial dependence on spins, i.e. we remove both 
the Jacobian J and the volume factors Vt- This means that the contribution from 6j'-symbols is effectively 1. We 
approximated the 9j-symbols by 6j'-symbols, provided the triangle inequality (|33[) is satisfied. Thus, the 9j'-amplitudes 
become trivial as well. 

Before we write down this approximation, we change to a notation in terms of indices — writing (x/i) instead of e 
and (xfih / ) in place of /. The m e 's that were associated to edges of n* are designated by (x/i), i.e. we write m I(1 . The 
m e 's that were associated to diagonal edges are designated by the plaquette {xpLv) in which the diagonal edge lies, 



We will ignore additional factors that would come from integrating over the quadratic order in fluctuation variables. 
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and we write them as fh X n V . The final result is 



{m xp } {m xp<T } 



x exp 



ini \ V p Lp x + b xp \ m xp 



47ri \ V p ip x + b xp + V a (p x + b xa \fh xprj 



(52) 



In the last term of the exponent, the repeated indices p and a are only summed over the pairs p < a. This expression 
is very similar to what we found in paper I. There are only three differences: 

• Here, we have a second set of monopole variables rh which arise from the discreteness of spins on diagonal edges. 

• The norm of J is j, while it is j + 1/2 in paper I. 

• The overall normalization is not clearly determined by the present derivation. 

As in paper I, we can factor off the Coulomb interaction between the currents: for this purpose, we choose the 
particular solution as 



Then, after changing variables 



and using the identity 



we obtain 



'•Pi 



A^v^ 



Yx 



(\XjWc) = 



(53) 
(54) 

(55) 



x exp 



( -p fx A( Px + 4iri | V M (<p x + &- x V u b xv ) + b x 



Ani |V P ((p x + A l V p b 

+ V„ (lfi x + A _1 Va&«7) + 



— T a A -1 T a 



(56) 



VI. SUMMARY AND DISCUSSION 



In this paper, we have derived an approximative model for 3-dimensional SU(2) lattice gauge theory. Its degrees 
of freedom can be interpreted as a dual gluon and monopole-like field. We propose it as a generalization of the 
photon-monopole representation of Polyakov [l3], [l4| and Banks, Myerson and Kogut (l5j . 

The derivation starts from the exact spin foam representation of SU(2) lattice gauge theory and extends an earlier 
work of Diakonov and Petrov (20| . The line of argument and the result are similar as in paper I. The crucial step is 
a stationary phase approximation: it can be viewed as the zeroth order of a weak-coupling perturbation theory that 
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has two types of non-trivial field configurations as a background: 1. the monopole-like excitations that arise from 
the compactness of the gauge group, and 2. the particular solution of the Gauss constraint which carries information 
about the large-distance defect created by the current. 

The resulting model contains two interaction terms: firstly, a current-current potential that is essentially the 
tree-level Coulomb interaction one would get from a purely perturbative treatment. Secondly, a coupling between 
monopole-like excitations, dual gluons and current. This coupling is similar to the photon-monopole coupling of U(l), 
but nonlinear. 

At the technical level, the present derivation is more complicated than in the first paper, and we have to make 
several heuristic arguments. To make this more precise, we would need the asymptotics of 9j-symbols, or another, 
systematic way of analyzing the sum over spin foam amplitudes. 

In spite of its heuristic character, the derivation from the spin foam representation is particularly interesting: it is 
the representation in which the strong-coupling expansion is performed, and in this regime confinement is a direct 
consequence of the spin foam diagrams. By analyzing the same representation at weak coupling we can ask the 
question how this stringy behaviour is preserved or modified as we go to the continuum limit. 

As we translated the steps of paper I to the spin foam representation, we found more solutions to the stationary 
phase condition than in the first paper. We learned that the Gauss constraint of the form (|4"Tj) covers only a subset 
of solutions and cannot capture all aspects of the non-abelian theory. When the spin of the Wilson loop is integer, 
the larger set of solutions includes spin foams of the type that are responsible for color screening at strong coupling. 
We suspect them to be also the cause for color screening in the continuum limit. 

The model of this paper can be tested by simulations and compared to simulations of the full theory. By summing 
over the monopole variables in expression (|52p one can remove the phase factors and translate them into a constraint. 
The latter can be enforced by a Gaussian damping factor. 
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